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ABSTRACT 



Aims. Emergent model spectra of neutron star atmospheres are widely used to fit the observed soft X-ray spectra of different types of isolated 
neutron stars. We investigate the effect of Compton scattering on the emergent spectra of hot (Jeff ^ 10* K) isolated neutron stars with weak 
magnetic fields. 

Methods. In order to compute model atmospheres in hydrostatic and radiative equilibrium we solve the radiation transfer equation with the 
Kompaneets operator. We calculate a set of models with effective temperatures in the range 1 - 5 -10* K, with two values of surface gravity 
(log g = 13.9 and 14.3) and different chemical compositions. 

Results. Radiation spectra computed with Compton scattering are softer than those computed without Compton scattering at high energies 
(£ > 5 keV) for light elements (H or He) model atmospheres. The Compton effect is more significant in H model atmospheres and models with 
low surface gravity. The emergent spectra of the hottest (^eff > 3 • 10* K) model atmospheres can be described by diluted blackbody spectra 
with hardness factors ~ 1.6 - 1.9. Compton scattering is less important in models with solar abundance of heavy elements. 
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1. Introduction 

Relatively young neutron stars (NSs) with ages < 10* yr 
are hot enough (Tes ~ 10* K) and can be observed as 
soft X-ray sources. Indeed, the thermal radiation of the iso- 
lated NSs was discover ed by the X-ray observatories Einstein 
and EXOSAT {Chens. & Helfandl 1 9831: lBrintanann'& Qgelman 
1987HCordova et al.«1989lK At the present time, the thermal ra- 
diation of a few tens of isolated NSs of different kinds, from 
anomalous X-ray pulsars to millisecond pulsars, are detected. 
The thermal spectra of these objects can be described by black- 
body spectra with (color) temperatu res from 40 to 700 eV (see, 
for example. iMereghetti et al.l2002l) . 

The nature of isolated NS surface layers is not known 
exactly. Under some conditions (depending on surface tem- 
perature, magnetic field strength, and chemical composition), 
a surface can be solid, liquid, or have a plasma envelope 
(|Lai & Salpeter 1997; Lai 2001). In the last case, the enve- 
lope can be considered as an NS atmosphere, and the structure 
and emergent spectrum of this atmosphere can be computed 
by us ing stellar model atmosphere methods (e.g. Mihalas 
19781) . Such modeling has been performed by many scien- 



tific groups, beginning with iRomanil (119871) . fo r an isolated 



NS model atmospheres without a magnetic field (Zavlin 
19961 iRaiagopal & Romanil IT996t IWerner & Deetiea .2001 



e t alj 



Gansicke et al ■I I2OO2I: IPons et alJl2002hT as well as for models 
with strong (g > lO'^ G) magneti c fields (Shibanov et al.'l9^ 
Raiagopal et alJl996l:IOzeJ200ll:lHo & LaL200L 2003. .20041) . 



These model spectra were u sed to fit the observ ed isolated NS 
X-ray spectra (see review bv lPavlov et al ] l2002h . 
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One of the important results of these works is as follows. 
Emergent model spectra of light elements (hydrogen and he- 
lium) NS atmospheres with low magnetic field are signifi- 
cantly harder than the blackbody spectra of the temperature 
Teff. These elements are fully ionized in atmospheres with 
Teff > 10* K. Therefore, the true opacity in these atmospheres 
(mainly due to free-free transitions) decreases with photon en- 
ergy as E^^. At high energies electron scattering is larger than 
the true opacity and photons emitted deep in the atmosphere 
(where T > T^ff) escape after a few scatterings on electrons. 
In the all of previous works, the model spectra of isolated NS 
were calculated with coherent (Thomson) electron scattering 
taken into consideration. As a result, hard photons, which are 
emitted in the deep hot layers of the atmosphere, escape with- 
out changing their energy. But if we take Compton scattering 
into account, the hard photons will lose energy at each scatter- 
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ing event. Therefore, such Compton down-scattering can affect 
emergent spectra of Hght elements model atmospheres of iso- 
lated NS. 

It is well known that the Compton down-scattering deter- 
mines the shape of emergent model spectra of hotter NS at- 
mosp heres with Ts-a ~ 2 ■ 10^ K and close t o the Eddington 



Umit dLondon et all 1 9861: iLapidus et all 1 9861: [Ebisuzakil 1987 



ZavUn & Shibanovl 1991; Pavlov et al. 



19911) . These model 



spectra describe the observed X-ray spectra of X-ray bursting 
NS in low-mass X-ray binaries (LMXBs), and they are close 
to d iluted blackbo dy spectra with a hardness factor fc ~ 1.5 - 
1.9 (lLondon etal.lil986: Lapidusetal 1986: Ebisuzakill 19*87 



atmospheres with Compton scattering taken into account are 
not calculated for relatively cool atmospheres with T^ff < 10^ 
K. Therefore, at present, the effect of Compton scattering on 
the emergent spectra of isolated NS model atmospheres with 
Jeff < 5 ■ 10^ K is not well known. It is necessary to point out, 
that the diluted blackbody spectrum has the same spectral en- 
ergy distribution as the black body spectrum with a given tem- 
perature Tc, but it has a lower flux. As a result, the bolometric 
flux of the diluted blackbody spectrum is lower than the bolo- 
metric flux of the blackbody spectrum with temperature . The 
bolometric flux of the diluted blackbody spectrum coiTesponds 
to the effective temperature T^ff < T^. The hardness factor is 
determined as the ratio of these temperatures: /c = T^jT^ff. 

In this paper, we compute model atmospheres of NSs with 
Compton scattering taken into consideration and investigate the 
Compton effect on the emergent model spectra of these atmo- 
spheres. We consider the importance of Compton scattering 
qualitatively in ^ Our methods of calculation are outlined in 
^ while results and conclusions are briefly discussed in §^ 
andm 

2. Importance of Compton scattering 

First of all, we consider the Compton scattering effect on emer- 
gent model spectra of isolated NS atmospheres qualitatively. 
It is well known that in the non-relativistic approximation 
{hv, kTs « nigC^) the relative photon energy lost due to a scat- 
tering event on a cool electron is 



^E 



hv 



nieC 



.2 ■ 



(1) 



Each scattering event changes the relative photon energy by 
this value. It is clear that the Compton scattering effect can 
be significant, if the final photon energy change is compara- 
ble to the initial photon energy. Therefore, we can define the 
Comp tonization parameter Zcomp (see also ISuleimanov et al.l 
2006 



Zc, 



hv 



omp 



max((T*)^T,*), 



(2) 



where max((T*)^, t*) is the number of scattering events that the 
photon undergoes before escaping, and t* is the Thomson op- 
tical depth, corresponding to the depth where escaping pho- 
tons of a given frequency are emitted. We can expect that 
Compton effects on emergent spectra of NS model atmospheres 
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Fig. 1. Comptonization parameter Zcomp vs. photon energy for 
neutron star model atmospheres with different chemical com- 
positions. 

are si gnificant if the Comptoniz ation parameter approaches 
unity ( Rvbicki & LightmanI 19791) . Because of this we compute 
Zcomp at different photon energies (see Fig.[T]| for hot NS model 
atmospheres with different chemical compositions. These mod- 
els were computed by using the method described in the next 
section, with Thomson electron scattering. It is seen from Fig.[T] 
that the Comptonization parameter is larger (0.1 - 1) at high 
photon energies (E > 4-5 keV) for H and He model at- 
mospheres. Therefore, we can expect a significant effect of 
Compton scattering on the emergent spectra of these models. 
On the other hand, Zcomp is small for the model with solar 
chemical composition of heavy elements. The Compton scat- 
tering effect on the emergent spectrum of this model has to be 
weak. 

This qualitative analysis shows that Compton scattering can 
be significant for fight element model atmospheres of isolated 
NS and we investigated this quantitatively in more detail. 

3. Method of calculations 

We computed model atmospheres of hot, isolated NSs subject 
to the constraints of hydrostatic and radiative equilib rium as- 
sumin g planar geometry using standard methods (e.g. iMihalas 
1978h . 

The model atmosphere structure for an NS with effective 
temperature Teff , surface gravity g, and given chemical compo- 
sition is described by a set of differential equations. The first is 
the hydrostatic equilibrium equation 



—— ^g-An \ Hy dv, 

dm Jo c 



where 



GMns 



4s VI -Rs/Rns' 



(3) 



(4) 



and when Rs = IGMf^^/c^ is the Schwarzschild radius of the 
NS, ky opacity per unit mass due to free-free, bound-free, and 
bound-bound transitions, cr^ the electron (Thomson) opacity. 



Suleimanov & Werner: Compton scattering in isolated neutron stars 



3 



Hy Eddington flux, f g a gas pressure, and the column density 
m is determined as 



dm 



-pdz. 



(5) 



The variable p denotes the gas density and z is the vertical dis- 
tance. 

The second is the radiation transfer equation with 
Compton s cattering take n into a cc ount using the Korn p aneet s 
operator ("Komoaneets' 'l 957t IZavhn & Shibanovl Il99ll: 
[Orebenev & Sunyaev 2002): 



(6) 



d I dJy Teff , / CJyS\ 



where x - hvlkT^s is the dimensionless frequency, fv{Ty) the 
variable Eddington factor, Jy the mean intensity of radiation. By 
the black body (Planck) intensity, T the local electron temper- 
ature, and C - c^h^ / 2(kTeff)^. The optical depth Ty is defined 

as 



dzy - (ky -H CTg) dm. 



(7) 



These equations have to be completed both by the energy bal- 
ance equation 



Jo 



ky (Jy - By) dV " CTg 



kT 



(8) 



^4^^ JydV-^j^^ Xjy(l + ^)dv'^^Q, 

the ideal gas law 



^ N,„,kT, 



(9) 



where A^tot is the number density of all particles, and also by 
the particle and charge conservation equations. We assume lo- 
cal thermodynamical equilibrium (LTE) in our calculations, so 
the number densities of all ionization and excitation states of 
all elements were calculated using Boltzmann and Saha equa- 
tions. We take the pressure ionization efifects into account on 
atomic populations using the occupation probability formal- 
ism (H ummer & Mihalas 1988) as described by Hubenv et al. 
(Il994l) . 

For solving the above equations and computing the model 
atmosphere, w e used a version of the computer code ATLAS 
jKurucz 1 19701 Il993h . modified to deal with high tempera- 
tures. In particular, we take free-free and bound-free transi- 
tions into consideration for all ions of the 15 most abundant 
chemical elements, an d about 25 000 spectral lines of these ions 
are als o included; see 'ibrag imov et al] ( 2003 ) and Swartz et al 



( 120021) for further details. This code was also modified to ac- 
count for Compton sca ttering ( Suleimanov & PoutanenI 2006t 
Suleimanov et al ]l2006l) . 

The scheme of calculations is as follows. First of all, the 
input parameters of the model atmosphere are defined: the ef- 
fective temperature Teff, surface gravity g, and the chemical 
composition. Then a starting model using a grey temperature 
distribution is calculated. The calculations are performed with 



a set of 98 depth points mi distributed logarithmically in equal 
steps from m ^ 10"^ g cm"^ to mmax- The appropriate value of 
found from the condition ^J Ty^b-i,f-f(mmax)Ty(min^^) > 
1 at all frequencies. Here Tyb f.f f is the optical depth com- 
puted with only the true opacity (bound-free and free-free tran- 
sitions, without scattering) taken into consideration. Satisfying 
this equation is necessary for the inner boundary condition of 
the radiation transfer 

For the starting model, all number densities and opaci- 
ties at all depth points and all frequencies are calculated. We 
use 300 logarithmically equidistant frequency points in our 
computations. The radiation transfer equation (|6]l is n on-linear 
and is solved i terative ly by the Feautrier method (Mihalas 



1978, ; see also Zavlin & Shibanov 1991; Pavlov et al. 1991 



Grebenev & Sunyaev 2002). We use the last term of Eq. (|6]l 
in the form xJ'y(l + CJ'y~^ /x^), where 7J, ' is the mean inten- 
sity from the previous iteration. During the first iteration we 
take 7^"' = 0. Between iterations we calculate the variable 
Eddington factors fy and hy, using the formal solution of the 
radiation transfer equation in three angles at each frequency. In 
the considered models 4-5 iterations are suflicient for achieving 
convergence. 

We used the usual condition at the outer boundary 



tlyjy, 

OTy 



(10) 



where hy is the surface variable Eddington factor. The inner 
boundary condition is 



dJy 

dTy 



dBy 

dTy 



(11) 



The outer boundary condition is found from the lack of incom- 
ing radiation at the NS surface, and the inner boundary condi- 
tion is obtained from the diffusion approximation Jy x By and 

Hy ^ I 13 X dBy I dTy. 

The boundary conditions along the frequency axis are 



Jy — By 

at the lower frequency boundary (v^ 
kTf.f[), and 



(12) 



lO''* Hz, /!v™„ <K 



dJy . , 7 eft , l, CJy\ 



CJy 



at the upper frequency boundary (v„ 



(13) 



10'^ Hz, hVr, 



kTgff). Condition ( fT2l i means that at the lowest energies the true 
opacity dominates over scattering ky '» os, and therefore Jy ^ 
By. Condition (fT3T l means that there is no photon flux along the 
frequency axis at the highest energy. 

The solution of the radiative transfer equation (|6]l was 
checked for the energy balance equation ([8]), together with the 
surface flux condition 



47r 



Hy(m = 0)dv = crT^g = 47r//o. 



The relative flux error 
Ho 



SH(m) = 1 - 



Q Hy(m)dv 



(14) 



(15) 
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and the energy balance eiTor as functions of depth 

kT 



Jo 



ky (Jy - By) " (T g J X 



(16) 



' eff 



14 j JydV-'-y- J xJy{\+^)dV 



were calculated, where Hy(m) is radiation flux at any given 
depth m. This quantity is found from the first moment of the 
radiation transfer equation: 



dfyJy 

dTy 



= Hy 



(17) 



Temperature corrections were then evaluated using three dif- 
ferent procedures. The first is the integral A-iteration method, 
modified for Compton scattering, based on the energy balance 
equation (O. In this method the temperature correction for a 
particular depth is found from 



X°° [(^v diag - 1)/(1 - CyAv diag)] ky{dBy/dT) dv 



(18) 



Here Uy - (Tsliky + cr^), and hy diag is the diagonal matrix 
element of the A operator (see details in Kurucz 1970l) . This 
procedure is used in the upper atmospheric layers. The second 
procedure is the Avrett-Krook flux correction, which uses the 
relative flux error e/f(m) and is performed in the deep layers. 
And the third one is the surface correction, which is based on 



the emergent flux error. See iKuruca (1 197ft) for a detailed de- 
scription of the methods. 

The iteration procedure is repeated until the relative flux 
error is smaller than 1 %, and the relative flux derivative error is 
smaller than 0.01 %. As a result of these calculations, we obtain 
a self-consistent isolated NS model atmosphere, together with 
the emergent spectrum of radiation. 

Our method of calculation was tested on a model for X- 



ray b ursting NS atmospheres dPavlov et alj|1991l; iMadei et al 
I2004i) . We found that our models are in good agreement with 
these calculations. 



4. Results 

We use this method to calculate a set of hydrogen and helium 
isolated NS model atmospheres with effective temperatures 1, 
2, 3, 5 -10* K and surface gravities log g - 13.9 and 14.3 
was calculated. Models with Compton scattering and Thomson 
scattering were computed for comparison. Part of the results 
are presented in Figs.|2]-|4] 

The Compton effect is significant for the spectra of hot 
(Teff > 3 ■ 10^ K) hydrogen model atmospheres at high ener- 
gies (Fig.lU. The hard emergent photons lose energy and heat 
the upper layers of the atmospheres due to interactions with 
electrons. As a result, the high-energy tails of the emergent 
spectra become similar to Wien spectra. The temperature also 
increases in the upper layers of the model atmospheres and 
chromosphere-like structures appear. Temperatures of these 
structures are close to the color temperatures of the Wien tails 
of the emergent spectra. Moreover, the overall emergent model 
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Fig. 2. Top panel: Emergent (unredshifted) spectra of pure 
hydrogen, low gravity (log ^=13.9) NS model atmospheres. 
Solid curves - with Compton effect; dashed curves - without 
Compton effect; dotted curves - blackbody spectra; thin dash- 
dotted curves - diluted blackbody spectra with hardness fac- 
tors 2.79, 1.87, and 1.6 for models with Teff = 1,3 and S-IO*" 
K. Bottom panel: Temperature structures of the correspond- 
ing model atmospheres. Effective and color temperatures are 
shown by dotted and dash-dotted lines, respectively. 

spectra of high temperature atmospheres in a first approxima- 
tion can be presented as diluted blackbody spectra with color 
temperatures that are close to Wien tail color temperatures 



Jc 



fcT, 



eff, 



(19) 



where /c is hardness factor and is connected to Hy as fol- 
lows: 



dv dv 
4KHy — , — 
dE dE 



1.6022- 10-^erg/keV 
h 



These results are similar to those obtained for model at 
mosphere s and emergent sp ectra of X-ray bursti ng NS in 
LMXBs dLondon et al] Il986l; [Lapidus etall Il986l; lEbisuzaki 



19871; iMadeill 1 99 ll; IPavlov eTakll 1 99 ll iMadei et al.l2004 

The Compton scattering effect on the emergent model spec- 
tra of high gravity atmospheres is less significant (Fig.O. The 
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Table 1. The hardness factors fc for hydrogen and helium NS model atmospheres with different log g computed for models with 
Compton (C) and Thomson (T) scattering. The values in the brackets show the deviations e of the model spectra from diluted 
blackbody spectra with color temperature Tc = fcT^fs (see Eq.l20]l. 
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Fig. 3. Top panel: Emergent (unredshifted) spectra of pure 
hydrogen NS model atmospheres with different surface grav- 
ities (solid curves - high gravity (\ogg - 14.3) models; dashed 
curves - low gravity models (log^ = 13.9)). For comparison the 
model spectra without Compton effect are shown for the hottest 
high gravity model (dotted curve). Bottom panel: Temperature 
structures of the corresponding model atmospheres. 



reason is a relatively small contribution of electron scattering to 
the total opacity in high gravity atmospheres compared to low 



gravity ones. The mass density in the high gravity models is 
higher, and the opacity coefficient (in cm^/g) is independent of 
the density for electron scattering and proportional to the den- 
sity for free-free transitions. In the presented models, H and 
He are practically fully ionized. Therefore, free-free transitions 
dominate the true opacity. 

The Compton scattering effect on helium model atmo- 
spheres is also less significant than on hydrogen model atmo- 
spheres with the same Teff and log g (Fig.|4l). The reason is the 
same as in the case of high gravity models. The ratio of electron 
scattering to a true opacity is smaller in the helium models. It 
is interesting to notice that, in the case of Thomson scattering, 
the emergent model spectra of helium atmospheres are softer 
than those of the hydrogen atmospheres. But in models where 
Compton scattering is taken into consideration, the spectra of 
helium atmospheres are harder than the spectra of hydrogen 
models. 

We fitted the calculated emergent spectra of this set of mod- 
els by the diluted blackbody spectra in energy range 0.2 - 10 
keV. We computed the value 

v2 



1 - 



^E,j(/c7"eff) 



(20) 



as a measure for the deviation of the model spectrum from 
the diluted blackbody spectrum B^. Here is the number 
of energy points E,j between 0.2 and 10 keV. CoiTesponding 
hardness factors are presented in the Table 1 together with e. 
Spectra of the hot models with Compton scattering are better 
approximated by a diluted black body, especially for hydrogen 
models. But the hardness factors are changing only slightly. 
Data from Table 1 can be used for interpretating of observed 
isolated NS X-ray spectra. The hardness factor can be found 
from these data if the observed color temperature T"^ and grav- 
itational redshift z of the NS are known. As soon as /c is de- 
termined, the apparent NS radius Roc can be found from the 
observed flux /e with the relation: 



(21) 



where d is distance to the NS. 

We also computed one isolated NS model atmosphere with 
solar chemical abundance and T^f^ - 3 ■ 10^ K and log §=14.3 
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Fig. 4. Top panel: Emergent (unredshifted) model spectra 
of pure helium, low gravity (logg - 13.9) NS model atmo- 
spheres. For comparison the model spectra of pure hydro- 
gen atmospheres are shown by dashed curves. Model spec- 
tra of the hottest atmospheres without Compton effect are 
shown by dash-dotted and dash-dot-dotted curves. Bottom 
panel: Temperature structures of the corresponding model at- 
mospheres. 

(see Fig.|5]l. The model was calculated with Thomson and 
Compton scattering and we found that the Compton effect on 
the emergent spectrum is very small. 

In our calculations we ignored deviations from the ideal 
gas EOS dU other than the ionization by pressure. We checked 
this approximation as follows. We found gas pressures from the 
temperatures and densities at all depth points for a pure helium 
atmosphere with Teff = 10^ K and logg = 14.3 by interpola- 
tion from the OPAL EOS tablefl We compared the obtained 
gas pressures with the values calculated in our model. The ra- 
tio of these two gas pressures are shown in Fig.|6] We find that 
the OPAL EOS gas pressure is equal to the gas pressure in our 
model with an accuracy of about 1 % up to column density ~ 
10^ g cm"^. At larger depths the OPAL EOS gas pressure is sig- 
nificantly larger due to partial electron degeneracy. All of the 
emergent photons are emitted at column densities lower than 

' http://www-phys.llnl.gov/Research/OPAL/ 
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Fig. 5. Emergent (unredshifted) model spectra of high gravity 
(logg - 14.3) NS atmospheres with the solar abundance of 15 
most abundant heavy elements with (dashed curve) and with- 
out (solid curve) Compton scattering. The dotted curve is the 
corresponding blackbody spectrum. 




10' lO'^ 10^ 10' 10' 10^ 10° io' 
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Fig. 6. Ratio of the gas pressure obtained by interpolation from 
OPAL EOS tables (see text) to the gas pressure from our cal- 
culations for the pure He model with Teif = 10^ K and logg - 
14.3. 



10^ g cm"^ (for example, photons with energy 10 keV are emit- 
ted at column density ~ 20 g cm"^). The distinction between 
the OPAL EOS and (|9]) must be smaller for models with larger 
Teff and for pure hydrogen models. Therefore, the EOS devia- 
tion from (|9]l at the deepest atmospheric layers cannot change 
the emergent model spectra that we calculated here. 

5. Conclusions 

We have presented the results of calculations for the hot model 
atmospheres of isolated NSs with low magnetic fields and dif- 
ferent chemical compositions. The Compton effect is taken into 
account. We investigated the importance of Compton scattering 
for the emergent spectra of these models. The main conclusions 
follow. 
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Emergent model spectra of hydrogen and helium NS atmo- 
spheres with Teif > 1 • 10^ K are changed by the Compton 
effect at high energies (E > 5 keV), and spectra of the hottest 
(Teff > 3 ■ 10*" K) model atmospheres can be described by di- 
luted blackbody spectra with hardness factors ~ 1.6 - 1.9. At 
the same time, however, the spectral energy distribution (SED) 
of these models are not significantly changed at the maximum 
of the SED (at energies 1-3 keV), and effects on the color tem- 
peratures are not strong. 

The Compton effect is the most significant for hydrogen 
model atmospheres and in low gravity models. Emergent model 
spectra of NS atmospheres with solar metal abundances are af- 
fected by Compton effects only very slightly. 

Acknowledgements. VS thanks DFG for financial support (grant We 
1312/35-1) and the Russian FBR (grant 05-02-17744) for partial sup- 
port of this investigation. 
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